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We make a systematic theoretical analysis of the masses of the fermions and their variation with 
energy by solving the one loop renormalization group equation (RGB) in the Minimal Supersym- 
, metric Standard Model (MSSM). A deceptibly simple common mass for all the fermions around 115 

' GeV at GUT scale has been found by Deo and Maharana. Here we undertake the unfinished but 

' the important task of calculating the electroweak masses of the fermions at different energies. The 

OA , proposed parametric unification mass and group theoretic constants for the model are known. The 

' mass of the top quark and its descent is studied by an approximate method very carefully. We find 

' that the Ramond, Roberts and Ross value of the Wolfenstein parameter is reproduced and is nearly 

equal to 0.22. When raised to integral powers and multiplied by 115 GeV, the whole mass spectra 
of the remaining eleven fermions are obtained within experimental errors. We deduce formulae for 
the masses and plot them for all the 12 fermions from t — log-^ = to tx = 33; the GUT mass 
being Mx = 2.2 x 10^'' GeV. 
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ir^ , 1. INTRODUCTION 

O . 

• Pati and SalamTl pioneered the idea that leptons are the fourth colour; quarks and leptons should be brought 
Oh, under the same umbrella of one group so that all forces ( except gravity) can be understood in terms of one unifying 
force parameter near the Planck scale. The six leptons should be treated at par with the six coloured quarks. 
^ ' The elementary constituents of matter would become twelve only. With the availability of enormous data from 
, the high energy accelerators, phenomenological analysis backed by imaginative theories, it has been found that the 
' familiar standard model, which is the product group S'[/c'(3) S'J7l(2) (g) [/^(l), can explain all of them successfully. 
. ^ , This standard model is characterised by three coupling strengths of the weak, the electromagnetic and the strong 
■ interactions. It was conceived that all the three couplings should run, i.e., they change with energy, and eventually 
H ' become one at the grand unified scale. This could be made possible by solving the relevant RGE. The coefficients of 
. . 1 the theory are specified by the group structure alone. Many theoretical models were investigated and only a few years 
back, it has become essential to bring in supersymmetry. In the minimal version of the MSSM, the beta coefficients 
are such that they run the coupfing strengths to a single value ^ = acuT — ^ at a mass of Mx = 2.2 x 10^^ GeV. 

This has been the most attractive result in the current investigations in gauge theories 0- 

Next in order, the major challenge in particle physics, was and is the theoretical derivation of the mass spectrum 
of the quarks and leptons in the same successful MSSM theory. In this model, all the masses of the fermions and the 
mixing angles were being chosen arbitrarily to account for the 19 free parameters of the theory. In the absence of 
such a fundamental theory, it has been in vogue to pursue the method which has been known as 'textures analysis'. 
After finding a suitable texture, one can investigate further and possibly obtain unification mass parameters for all 
the fermions as a generalisation of the hypothesis by Georgi et alQ- Eventually, one can predict their individual 
masses by using the analytical MSSM group coefficients for a RGE for the masses of the fermions. As yet, this has 
not been successful. 

One of the 13 parameters of the standard model was first predicted in 1974 by Gaillard and LeeQ. Then came 
the work of the popular mass matrix ansatz of Fritzsch^^. A complete listing of textures and their relevance to 
experimental findings were made by Ramond, Roberts and Ross ( RRR)0. Since they also took the help of RGE, 
we shall present the essential results of their work in this section. Before that, we present the MSSM Lagrangian 
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QlMu^uUr + QLMo^dDR + IlMe^clEr + ulMn^uVr + h.c. (1.1) 
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We ignore the sparticles and consider the rigid part of the Lagrangian given by Demir • 

In the renormahsation schemes, the Yukawa couplings Mp{t) and v.e.v. vp{t) change with energy separately. The 
Dirac masses are given by 

mpit) = VF{t)MF{t) (1.2) 

The masses considered above are not the masses of the flavor eigenstates of the model. 
The one loop RGE, as written by Grzadkowski, Lindner and Theisen @ for Myit) are 

l6n^^^^i-GY+TY + SY)MY{t), (1.3) 

t — logln/Mz), /i is the renormahsation point and AIz is the mass of the Z-boson. Here Y=U stands for U-quarks 
(F=l,2,3): up(u), charm(c) and top(t); Y=:D stands for Down quarks (F=4,5,6): bottom(b), strange(s) and down(d); 
Y=N stands for Neutrinos (F=7,8,9): v^, f^, and i^r, and Y=E stands for electrons(F=10,ll,12): e, /x and r. G'y(i) 
contains the gauge coupling terms, given in Section-2 and 

Tu ^Tr(iMuMl+MNMl,y To ^ Tr [iMoM], + MeMI^ Tu =Tn,Te = Td, (1.4) 

and 

Su ^3MuMI+MdMI, Sd =3MdMI + MuMI, Se ^ SMeMI + MnM^, Sn ^ ^MnM]^ + MeMI. [l.b) 

To find the couplings, we have to solve the twelve differential equations and determine all the couplings/masses 
from only one value of coupling at t = tx or Mf{Mx)- First, we turn our attention to the mass of the top quark. 
Incidentally, Pendleton and Ross 0, Faraggifioj. and some others have predicted the value of the mass of the top 
quark which was around 175 GeV, even before the top was discovered. 

In this paper, all the one loop equations are solved following the same method used for gauge coupling RGE. A 
particular Mf{Mz) is obtained relating to Mf{Mx) and other fermions. Considering the top, first we make 'heavy 
top integral' approximation stated below. Eventhough we aim at single input value Mf{Mx) — Mu , which is 
independent of F, the integrals of the solutions for the masses are such that 

J Mlp{T)dT » j Ml^^^^{T)dT »> j Ml^^^,^{T)dT (1.6) 

Neglect of the integrals other than the top quark, gives a very good result for the Mtop{Mx) — Mu- Following this 
approximation procedure, we find from other eleven equations that Mtop{Mx) — Mf{Mx) — Mu — 115 GeV for all 
the fermions. 

We choose Mu — 115 GeV as the only input. Furthermore, MSSM has two Higgs. 

<(j)u>^Vu{t) = v{t)Sinfi{t) (1.7) 
< Vd{t) - v{t)CosP{t) (1.8) 

v\t) = vl{t)+vl{t) (1.9) 

tanm = (1.10) 

It is to be noted that v{t) is the vacuum expectation value of single Higgs of the Standard Model. v{Mz)=24:6 GeV, 
SinPsM{Mz) = We note that = ^ GeV=174 GeV, which is close to MtopiMz) = 175 GeV. For simplicity, 
we shall use this as the unit of energy whenever unspecified. 
The one loop RGE are 

levr^^^ = ( ^g! + Igl - TriZMuMl)] v^, (1.11) 



dt \ 20 



16-^^ = (^9l + l9l~Tri3MnMl,)\v,. (1.12) 



The calculation of 



(1.13) 
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in MSSM has attracted considerable attention, there exist extensive hterature, most of them arc given in reference p^. 
We note that below the mass scale Mz of the normal SM, one Higgs v.e.v. satisfies the one loop equation 



Wn^^^aiE. = (J^gj + _ Tr{iMuMl + iMoM],)^ vsm- 



(1.14) 



To maintain continuity, we assume that at t—0, 

tan(3{Mz) = tan(3sM{Mz) = 1. (1.15) 
Using the approximation given in equation H1.6|l . we get from equations Hl.ll|l . (|1.12l) and H1.13|l . 

tanm - exp (^-J^ [Mlpir) - MI,,^Jt)] dr^ . (1.16) 

Using the results given in scction-3, eauation l|3.27|) . we find that tan[3{t) is a slowly varying function of t and drops 
from 1 at 91 GeV to 0.9 at 10^^ GeV. So the calculation is much simplified if we take SinP = ^ = CosP. This is 
not ruled out by experiment!^. 

The authors (RRR), noted that there are only six possible forms of symmetric mass matrices with three non-zero 
eigenvalues and three texture zeroes, capable of describing the hierarchy of up or down quark mass matrices. Those 
are 

ag 

Ti = ( 6i I , T2 = ( a2 62 I , T3 = ( 63 

63 C3 

ae be 

T4 = ( 64 I , T5 = ( as 65 I and Te = ( ae I . (1.17) 

be C6 

Ts was the one first pioneered by FritzschjJ]. More recently Dimopoulos, Hall and RabyJJJ have analysed and included 
the leptons following GeorgiQ and solved the MSSM RGE with some degree of success. RRR tried to analyse all the 
cases and put them in a CKM matrix form, proposed by Wolfenstein 13] and diagonalise them to the texture types. 

C1C2 - siS2e"*"^ si + ciS2e""^ §2(53 - 54) \ 

VcKM = I -C1S2 - sie"*"^ -S1S2 + (C1C2C3C4 + siS4)e"'"^ S3 - S4 (1-18) 

Sl(s3-S4) -Cl(s3-S4) (C3C4 + SlS4)e*'^ / 

AV2 A A\^p + ir])\ 

-X I-AV2 AX^ . (1.19) 

AX^{1- p + ir]) -AX'^ 1 J 

Here Sj, q (i=l,..,4)are the sines and cosines of mixing angles. From the identities given by Dimopoulos, Hall and 
RabV[l^. following from equations H1.18|l and (|1.19|) 

X^ {sl + sl + S1S2 Cos4>y^^. (1.20) 

(f) is the CKM phase angle and 

MdY , fMuY ,2 fMdM. 

) = ^' «2 = ( — ) = A ; S4=[—^ 
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As will be discussed later, the Coscj) defined in we shall get 

Cos(/) = -A/2, (1.22) 

and 

S3~S4^ X^A{t) (1.23) 

where A(t) depends on t. 
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The small expansion parameter is A ~ 0.2 and A ~ 0.9 ± 0.1. Olcchowski and Poroski[^ were the first to write 
down the RG Equations for the parameters 

levr^^ = -3c{hj + hl)\J,,\, (1.24) 
dX 



0, (1.26) 

0, (1.27) 

0. (1.28) 

Here c = 2/3 for MSSM and Jcp is the irreducible phase of the CKM matrix. To arrive at equation (|1.25|) . we equate 

-1/7 



dt 
dp 

'dt 

drj 

'dt 



Ait) 



Mbottom{t)Mtop{t) 

MUMx) 



(1.29) 



top\ 

Using one loop RG equations H1.3|) to H1.5|l . we get 

A{Mx) = 1.00 and A{Mz) = 1.474. (1.30) 

RRR have made a complete listing and analysis. They arrived at the value A —==0.22 and the oft quoted result, 

rrir ■ rrip : rrie — mt, : nis : = 1 : A^ ; A'* 

rrit : nic '■ rriu = 1 : A'' : A^, (1-31) 

which is very well satisfied by experimentally found masses. This has been thought to be like a miracle. 

We shall take the masses of the twelve fermions to be real, wherever necessary. In section-2, we shall write the 
RG Equation with the MSSM coefficients and give the one loop exact solutions. In section-3, the unification mass 
for all fermions at GUT scale ranging from 113 GeV to 125 GeV as computed by Deo and Maharanajl^, will be 
discussed as a follow up of their letter. An expression for Wolfenstein parameter in terms of RGE coefficients is given 
in section-4. The reason to raise this parameter by integers np to obtain fermion masses except the top, Mp equaling 
A"^ multiplied by the parametric unification mass 115 GeV, is given. In sections-5,6 and 7, we suggest an alternative 
method of calculation of experimental masses of all fermions by a suitable self contained procedure which includes 
the gauge couplings as well. The equation for the running of all the masses of all the fermions is given in section-8. 
The results are given in tables and are also shown graphically for greater clarity. The concluding remarks are given 
in the section-9. 



2. RENORMALISATION GROUP EQUATIONS AND SOLUTIONS 



As stated, the gauge sector of the standard model is characterised by three coupling constants 33, 32 and gi of 
SUci'i) ® SUl{2) ® {/^(l), respectively. However, these couplings are not constants, they change with energy/mass 
values. The nature of variation is given by the solutions of the RG Equations. The coefficients are calculated by the 
specific nature of the Standard Model group. For MSSM, the three couplings at mass Mx = 2.2 x IOI6 GeVjil] unite 
to a unified coupling constant g^/47r=l/24.6. Here, the supersymmetry descends from Mx down to AIz — 91 GeV, 
as suggested by Witten. The RG Equations for the couplings in the lowest order are given by 

16n'^ = c^aUt), ^ = 1,2,3. (2.1) 

The coefficients are ci = 6.6, C2 =1, C3 = -3. The first two coefficients are positive, indicating that C/y(1) and 
SUl{2) are not asymptotically free, whereas the SUc{3) colour group is free and makes the entire product group 
asymptotically free . This implies that in a perturbative formulation, the higher order contributions are small and 
can be neglected. Therefore, we shall use eauation H2.1|l only in the gauge sector, with two HiggsfTsj. 
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The running parameter t is defined as t = loge-rj- so that it varies from to loge^rf^ — 33. The solution to RG 



Eauation (|2.1|l is 



47r 



All Ci 



91 



2lT 



(2.2) 



Here, g. 
and '"^ 



(7^(0) are the couphng strengths in the electroweak scale Mz- Taking the value of Mx = 2.2 x 10"'^^ GeV 
24.6, we calculate the values of ^ = 59.24, ^ = 29.85 and ^ = 8.85. These are consistent with the 

9p ' 91 ' 92 93 ^ 

experimental results. Thus the three coupling strengths are descendants of one coupling constant gy. 

Taking the clue from eauation H2.1|l . we rewrite the Yukawa sector SUSY RG Equations given earlier in equation 
(|1.3() . which had also been written by Babu^^ following Georgi and Glashow, and Eitchen et al^3| in the following 
way; 



167r 



di 



= ApMUt) + [Ypit) - GF{t)]MF{t) 

= AFMUt) + ZF{t)MF{t). 



(2.3) 
(2.4) 



The masses are in units of 175 GeV. We repeat for ready reference that the 12 fermions are suffixed as F = 1,2,... 12. F 
= 1,2,3 are the [/-quarks: top(i), charm(c) and up(w). Similarly, F=4,5,6 denote the D-quarks: bottom(6), strange(s) 
and down(d); F=7,8,9 are the i?-leptons: electron(e), muon(/i) and tau(T); and i^=10,ll,12 are the A^-neutrinoes 

V^, l^fj,, Vr- Further, Mi = Mtop, M2 = Mcharm, ■■■ , M12 = Mtau- 

Af is a group theoretic factor whose value is '6' for quarks, i.e., F=l,2,- • • ,6 and '4' for the leptons i.e., F=7,8,- • ■ 
, 12. The positive values indicate the field theory containing Yukawa couplings only and may not be asymptotically 
free. 

Yp is the mixing term which can be put in matrix form 



Yf = Y1 ^FHA4(t)Afff (t), = 1, 2, • • • , 12. 



(2.5) 



H 



In MSSM, the matrix Afh is specified by the 144 elements given below. 
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0/ 



(2.6) 



The diagonal elements of Afh have been taken out as the cubic term in equation |^J|; so they are zero. 

As the model is minimal supersynimetric, the gauge factors Gi?(i), which are the sum of gauge couplings, are fixed, 
we take the values from reference and |1^. 



13 



Guit) = ^glit) + iglit) + ^^gl{t) ^Y^^ug^t) 



15 



(2.7) 



7 16 

TT5?(i) + 3g2'W + ^53W, F = 4,5,6 



F =1,2,3 stand for U and they are degenerate electromagnetic gaugewise. Similarly, 

GE{t) = ^.g?(t) + 3<7i(t), F = 7,8,9 

5 

and GN{t) = -gl{t) + 3gl(t), F= 10, 11,12. 



(2.8) 
(2.9) 
(2.10) 
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Here = = 0, as the leptons do not have the strong colour interaction. We shall need the integrals, 



and 



3 



-i,/ =E^log(l-^^) (2.12) 



2 — 1 

= E — log%- (2-13) 

Deo and MaharanaT^ made a very important observation that eauation (|2.4l) . which has to be solved for a given 
fermion, does not contain the coeficients of the same mass in the matrix A of equation 1)2. 6|l . This fact has been 
overlooked by all previous authors. As a result, the calculational dedails become erroneous and the values obtained 
are unreliable. The present approach gives a hope of a simple method of entangling the mass due to finding the 
solution of 12 difFerentential equations. As such, the terms Zp{t) can be exponiented away. We introduce a subsidiary 
mass mpit), through 

Mpit) = mpit) exp (^j^ ^ Zf{t) dr^ , (2.14) 
such that Mp{Mz) ~ mF{Mz) = tof(O). They satisfy the equation 



167r^ = ^Fexp — / Zf{t) dr dt. (2.15) 

They look, astonishingly, similar to the gauge sector one loop RG equation (|2.1)l and can be solved exactly. Integrating 
eauation (|2.15|l from Mz to Mx i.e., from t^O to tx, we get 

Stt^ Stt^ f*^ / 1 /■* \ 

2,,r , = ^TTFT + / exp / ZF{T)dT . (2.16) 

Putting back the exponential, 

Ml,(M,) _ M-i,,Afe) / I ^ A, r" ^ / I r' ^ ^^^^^ 



Ml{Mz) MliMx) ySn^ Jq J Sn^ Jq \8t: 

This is the exact one loop solution. Mp{Mz) is the mass of the fermions at Mz- 

The descent or ascent running of the masses from GUT Mx to electroweak Mz, can be obtained by integrating 
eauation (|2.15|l from t = tx to t. The result which is not given in reference [l^] is 



(2.18) 



This is also one loop exact. By solving equations H2.17|l and H2.18(l . we can find Mf{Mx) and Mpit) repsectively. 



3. ORIGINAL MASS OF ALL FERMIONS AT Mx ? 



The gauge integrals over Gpit) are easily and acurately calculable. The most diSicult task is to evaluate YF{t). 
Even though, it does not contain Mpit), it is a sum of squares of moduli of the masses of all fermions. For example, 
from the matrix as given by eauation l|2.6(l . 

Yt,p{t) = SM^it) + 3M^{t) + Ml{t) + Ml (t) + Ml {t) + M^ (t). (3.1) 

There is mixing of six other fermions for the top. 
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In the 'top heavy integral' approximation, one retains only those terms containing M'^^p{t) occuring in any integral 
with Yp{t) . Consider the top case. Ytop{t) can be set equal to zero in the RG Equation for the top. The top mass is 
then given by a simple expression 



where, using integrals H2.11|l to H2.18|l . 



^*op = n % = 0.086 (3.3) 
,^1 \9uJ 

and (3.4) 

-a^i "'11(1-1;^) -"^""^ («) 

Putting these values, the original mass of the top, at an energy of 2.2 x 10^^ GeV was, approximately, 

Mtop{Mx) = Mtop{Mz){l- DtopT^'^Clil ^llA GeV. (3.6) 

This is Deo-Maharana result. 

We can write a general formula for all fermions 

MlpjMz) MlpiMz)^ 



where 



M'p{Mz) M'p{Mx) 



3 ^ 

It is not so easy to calculate Yp for other fermions even with 'heavy top integral' approximation given by equation 
(|1.6|) . So we consider the Yukawa-like coupling hp{t) for F=2,3, 12, which is related to Mp{t) as 



1 



Mpit) = hFit)exp ( / GF(r)dT ) . (3.10) 







The RG Equation for hp{t) is, 

Sir^d {log h%{t)) = AFhl{t)e^p(^~-^ Gp{T)dT^ dt + Yp{t)dt. (3.11) 

Here hp{t) contains the usual /?-angle factors of the two Higgs system. As has been discussed in Section-1, we may 
not need this angle in our approach to the problem at hand. If F > 1, the heaviest fermion next to the top is the 
bottom. The first term of eauation l|3.11(l . at the Mz scale, is 1/1235 times smaller, whereas the last term contains 
one htop in hbotm- To a good approximation, and to begin with, we shall neglect this term and express hpit) in terms 
ofYp{tys. Then 



^^^74(Mz) 



Yp{T)dT = STT^log J ' ■ (3.12) 



Here, h?p{Mz) = Mj,{Mz) by definition. In the miits of top mass 175 GeV, 

hUMx) = M|(Mx) exp (^-^ ^ Gp{T)dT^ = Ml^iMx) exp (^-^ ^ (T)dT^ (3.13) 



Mf(Mz) = 1. (3.14) 
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So, 



I YF{T)dT = -Stt^ log M|,(Afz) (3.15) 



Ml{Mz) 



exp(^ / Y,{r)dr] = (3.16) 



This is true as long as the first term of equation ((^^J is negligible. Since there is Mp in this term, it is much more 
justifiable to set Diepton — 0. The general equation H3.7|l for the leptons gives 

Ml,,AMz) - Afa^,„„(Af^)^'--"' ^'-^'^ 

where 



Cle.pt on — 



9l 



1=1 



^leptoni^x) — M^epton^^lepton (3.18) 

Putting the gauge constants and masses in the above, as experimentally reported, we get the masses at GUT scale 
for different leptons as given in Table-^ 



TABLE I: Unification scale mass for leptons 



Lepton 


Unification scale mass (GeV) 


e 


116 




116 


T 


116 




126 




126 




126 



At the grand unification mass, the electron, the muon and the tau climb to 116 GeV in this approximation. The 
calculation for the neutrinoes are not reliable as the isospin factors are uncertain and may be inaccurate. 

We are now left with the five quarks. For them, we attempt to find the next leading order approximation, i.e., we 
first set the first term equal to zero and obtain 

This is used in the calculation for Cq of eauation H3.8|l . which gives 



Using 



9 



„ _ 6 hlit) J^f 5fc.A - 



in eauation H3.9|l . we get 



87r2 7o M2 (Mz) l\ V 87r2 , 

We look for ways to calculate hqit). We can retain only the top quark in the RG equation in a slightly different way 
than what has been taken by RRrQ and get, 

87r2dlog/i^(t) ~ Nqdt. (3.24) 

We have 

^charm — 3, ^up — 3, -^bottom — f ; -^strange — ^down — 0- (3.25) 

To use this value in Yp, it is necessary to take an average as the couplings change very rapidly, 
or 

This expression for Yukawa coupling, is not much different from unity for all allowed t's. We arrive at the following 
result. 

The unification mass is calculated numerically from 

MgiMx) = Mtopl^^^] ' . (3.29) 



Q 



The results, for the quarks, are given in Table-ini We note that M^harm{Mx) = Mup[Mx) 7^ Mtop[Mx) and 

M,trange{Mx) = Mdown{Mx) ^ Mbottom{Mx). 



TABLE II: Unification scale mass for quarks 



Quark 


Unification scale mass (GeV) 


Top 


114 


Charm 


115 


Up 


115 


Bottom 


119 


Strange 


118 


Down 


118 



Thus, we have shown that all fermions seem to originate at an energy 2.2 x 10^^ GeV with equal mass of about 
115 GeV. Perhaps, this is due to the equality of A~ 1 in Wolfenstein's parametrisation of CKM matrix. In this 
perturbative method of solution for finding the unification mass, information about the masses of the 11 fermions has 
been lost due to cancellation of AP{Mz) in both r.h.s. and l.h.s. of the equation (|2.17() due to use of equations (|3.16() 
and H3.22|l . The result of a common mass at the origin is atleast a hypothesis and has an approximation |l4| . 

The top mass decreases with energy but all the other quarks and leptons, starting from the value at Mx acquire 
smaller and smaller values and become quite light in the electroweak scale. The descent or ascent eauation H2.18f) 
describing the 'run' can be put in a form like eauation (|3.7() . 

MlJMz) Mf„JMz) 
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where 



1 



Cpit) = apit) exp ( — / Yp{T)dT ] , (3.31) 



Dpit) - ^ ^ dhbpih) exp ^ ' Yp{T)dT ) , (3.32) 



g~cit 



(3.33) 



and hp(t) — JJ^ 



1 - 



(3.34) 



For the top, we can take — > and calculate variation of its mass from Mu ~ 115 GeV to the top mass 175 GeV 

Mtop{Mx) 



Mtopit) 



Ctop{t) + Dtop{t)j 

top / 



1/2 



(3.35) 



The values are given in Table-0of Scction-8. 

For other cases, we shall try to fit them into a scheme which is not only much simpler, but has better physical 
content. However, in the following section, we discuss a different route for solutions without specifying gauge factors 
completely. 



DEDUCTION OF WOLFENSTEIN PARAMETERS AND THE ROTATIONAL INTEGERS 



Let us construct a function Bp(t) such that 



Then 



Eauation (|2.17|l reduces to 



8TT^-logBUt) = Zpit). 



{sn^ Jo 



exp I -^^—^ / Zp{T)dT 



and exp Zp{T)dT 



BUMx) 
BUMzY 

BUMz) 
BUMx)' 

BUMz)' 

BUt) 
BUMxY 



Mp{Mz) 



1 



Ml(Mx) A. 



1/2- 



As indicated in the last section, except the top, we neglect the terms with Ap. Then 

If 



Mp{Mz) ~ Mp{Mx)^^Yrr\ = MF(Mx)exp 

Up[Mxj 



167r2 



= Mp{Mx)y' 



(4.1) 

(4.2) 
(4.3) 
(4.4) 
(4.5) 

(4.6) 



(4.7) 
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where the integral Ip is 

If = jj ZF{t)dt ^ jj ApGMl{t)MG{t) - Gpit)^ dt, (4.8) 

with AiG—0 to indicate that this is not for the top. Since the ratio of the two fermion masses can be put as -^j^ = A, 

the parameter M^, cancels out and A is the Wolfenstein parameter. The purpose of introducing the logarithms is to 
incorporate the possible multivaluedness of Ip and obtain the integral powers of A. From equation 1)4. 8|l . 

If ^ Zp{t)dt ^ r [Y,AFGMl{t)MG{t) - Gpit)] dt (4.9) 







G 



\jj{Y. ApGMlmiGit) - Gpit) + ApGMl{-t)MG{-t) - Gp{-t)^ dt 



(4.10) 



1 rtx dt I \ 

4 / dKT'^^^^^^^ [Y,ApGMUt)MG{t) - Gpit) + J2AFGMli~t)MGi^t) - Gpi-t) (4.11) 

J-tx ^ \ G G / 

[^"" dMpit) (Eg ApGMlit)MGit) - ^iGpit) + Gpi-t)) 



167 



(4.12) 



-Mu Mpit) {AFMlit)Mpit) -Gpit)^ 

where we have used equation (|2.4|l . Let us set 

Mpit) ^ Mc/e*^^^*^^^*""'^- = M[/e**^(*)"^, 
where np is an integer and we will call it rotational integer. Using the above in (|4.12l) . we get 

If = ^np^^ dQpit) ^^'^ ^ .. . (4.13) 

2 J-tx ^ ' Af + Y.gAfg-^ 

We have omitted inconsequential factor ^p- in the numerator. We note that for quarks Ap + Yl%=i~^~^'^~^^ ^'^^ 
for leptons 4+9=13. The integral is almost a constant and isolate the integer np, characterising F from the incgral. 
Using 1 = ^ XI i/ equation (|4.13|l , we have 

1 C_2 1 ptx 1 

^^=-^^I^EE^^W denit) ^. (4.14) 

^ H G -^-^^ Ah + }_^gAhg ~ mT 

In the above we have averaged over the twelve fermions. Retracing the steps and using Mnit) — Mije^'^'^^^^^'^u in 
eauation (|4.14|l . we finally get 

Ip^np^Y.Y.^HG j dt= ^nptxY^Y^^HG- (4.15) 

H G •'^ H G 

First, we shall be interested in the CKM matrix for quarks with lepton masses taken as zero. Then let G vary from 
1 to 6, whereas F will be taking values from 2 to 12, since all of them contains 6 quarks. Only the coefficients Ahg 
are needed to calculate Ip of equation (|4.15(l . The coefficients ApG have non-vanishing values for F=2,3,- • • ,12 and 
G=l,2,-.. ,6. 

For F = 2,---,6: Api + Ap2 + ■ ■ ■ + ApQ = 7 , (4.16) 

and 

Fori^ = 7,--- ,12: Api + Ap2 + ■ ■ ■ + ApQ ^ 9. (4.17) 

For the twelve fermions, the sum of the values of the coefficients Afg is (7 x 5) + (9 x 6) = 89. The average of Z, i.e., 
Z is geven 

89 

Z = 89/12, and Ip = rip tx —■ (4.18) 



12 



From cauation (|4.7|) . the Wolfenstein parameter A tmns out to be 

This is an excellent result in spite of the approximate estimates. The masses of all the fermions due to quark-lepton 
equivalence other than the top is 

Mf{Mz) ^ AfF(Mx)A"^ ~ y'Mu (4.20) 

The table 4IIll identifies the particles. We have increased np hy neighbourhood integers which we have called the 
rotational integers. 



TABLE III: Identification of fermions 





Mass in GeV 


Fermion 


2 


5.5 


bottom(b) 


3 


1.2 


chami(c),tau lepton(r) 


4 


0.264 


strange(s), muon(^) 


6 


0.012 


down(d),tau neutrino(!/T) 


7 


0.0028 


up(u) 


8 


6 X 10"'' 


electron(e) 


9 


1.33 X 10~* 


muon neutrino (j/p) 


16 


3 X 10"^ 


electron neutrino(i/e) 



5. GAUGE CONTRIBUTIONS 
5.1. The top mass 

It is easy to write down the equation for the top mass in terms of the mass Mp(AIx) and gauge couplings. 

1 rl \ 

^ ^ U-top 



Mtop = Mtop{Mx) [ ^ ) . (5.1) 

^top 



We shall also need the equation for the descent of top mass 



^^atop(i) + Aop(0- (5.2) 



The equation contains only the gauge factors. 

5.2. The masses of charm and up quarks. 

The RG Equation for the top is non-linear and intermixed. Specifically it is 

Sn^j^logMl^it) ^ -Gu + 6Mlj,{t) + Ytopit). (5.3) 

Following a texture analysis of the RG solutions, we introduce a function B{t),{ not to be confused with Bpit)), which 
satisfies 

STT^j^log B^{t) = -Gu + &Ml^{t). (5.4) 

Eauation H5.4|l does not specify the function B{t) completely except that d\ogB{t) = d\og{Mtop{t)) ■ This lone 
restriction gives us an infinite number of free choices. Because letting B— + ^B, where ^ is an arbitrary constant, does 
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not change the top mass function Mtop{t). For simpUcity, we take the second derivative of this function B{t) to be 
zero so that 



dlogB{t) = Cndt. 



(5.5) 



Cb is essentially d\ogMiop{t), Mtop{t) changes by only 30 to 50 GeV as the mass /i goes from 91 GcV to 2.2 x 10^^ 
GeV. Cb ^ (175 - (123 to 115))/175 ~ 0.297 ' to 0.342. We take Cb = 0.3. The nonequiness due to nonlinearity 
can also be seen as follows: 



dlogB{t) = dlogM{t) = dlog 



top 



dt 

Mo ^ 167r2 

Mo is an arbitrary constant. We can choose Mq suitably so that 



-Gu + Q 



Mi 



top 



Mo 



Mi 



top 



Mo 



Integrating from t = tx{Mx) to t = O(M^), we have 



1 



/167r2 ~ 0.3 = Ce 



-CBtx, 



and 



Bz 

B. 



u 



Furthermore, integrating from t = tx to arbitrary t, 

Bit) 



log 



Bu 



= (0.192)^ 



CB{t-tx), 



or 



m 

Bu 



= exp(Cs(t - tx)) 



Bz_ 
Bu 



1/6- 



Let us examine the case for the charm. We have 

d 



Sw'-logM^it) = -Gu + ^M^, 



= -Gu + l[6Mf-Gu] + lGu, 
= -iG. + 4.^|logB^(0. 



Integrating from tx to t, 



log 



M^HarJMx) 



ilog^-^ 
2 ^ Bl 167r2 



Gu{t)dt, 



tx 



Mcharmit) = Mu 
McharmiMz) 

and Mc 
From section-4, we now calculate 



m 

Bu 



1/2 



exp 



327r2 



Gu{T)dT 



tx 



Mui^y^Ci 



Bu 



1/6 



-1/4 

top 5 

3 



C,o'J^ = 1.8466, 



Ci 



-1/12 



top 



McharmiMx). 



A 



charm 



g-CBtx/6c>-l/12 ^ g_221 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



(5.10) 

(5.11) 
(5.12) 
(5.13) 
(5.14) 

(5.15) 
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Again this is a good result. The mass of the charm is 1.24 GeV from the equation 1)5.14(1 in very good agreement with 
the experimental value. The rotational integer np is three as in Table- HllI 

We continue further and consider the up quark. The 'heavy top integral' approximation of the RGE for the up 
quark is 



87r2- logM^ = -Gu + 3A4p(t). 



(5.16) 



Guessing from the general rotational integral parametrization, as shown in Table- HTTl as A^, we write this equation as 



Integrating, we get 



^Gu + -{6Ml,{t)), 



Mu 



Gjj 



El 



?7r2|logB2(t)_^8^2^1ogA4^(i). 



d 



\Mtop 
Bz\ f Mu 



1/2 



^1/4 



BuJ \Mt 



top 



1/2 



-1/4 



= 0.220. 



(5.17) 

(5.18) 
(5.19) 



This gives M^p — 0.0029 GeV, with A% = Mtop{Mx), a little lower value but quite close to the experimental 
result. The values depend crucially on Cb, which is itself not exact. It is also true that Mjj may be a little different 
from 115 GeV. 



6. MASSES OF THE DOWN QUARKS 

The calculation of the mass of up-quark has set the method of finding solutions close to the experimental values 
by using the texture analysis function B{t) ignoring the terms with coefficients Ap. We present a general recipe from 
the 'heavy top integral' approximation. Let the RG Equation for any F 7^ 1 be 



Stt^ ^ log Af| ^~Gp + NpMlp. (6. 1) 



Here Np can be zero as well. 



NpMl^it) = ^{6Ml^it)~Gu) + ^Gu, 



Np + 'i]p 
6 



8.^1 log i?^(i)) - (sn'^^logMl^it)^ + ^Gu. (6.2) 



r]p can be any arbitrary coefficient. All of them will satisfy the RG Equation in the 'heavy top integral' approximation. 
But they should let np and A be such that they are within rotational integers and gauge interaction contributions. 
Integrating from the known values, tx to 0, we get 



167r2 Jq " ' ' 967r2 " ' ' J \Bu J \Mi 



Mp{Mz) - Mt/exp ( / Gp{t) dt - / Gv{t) dt\\^\ ( --^ ) . (6.3) 
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For the down quarks, 



So 



8.^1 log Af.^ 

STT^^logMj 

8^2^1ogAf2 



-Gd + Mf^p, 

-Gd, 

-Gd, 



-Gd + 2M4^op - Mtop, 



-Gd + ^ml^ - Gu) - ^ml^ -Gu) + ^ 



Gjj 



6 6 

2.5 , 9 rf , 1.5 , n d , 9 , GjT 

-(8.^-logB^)--(8.^-logM^,) + ^. 



Integrating from ijf to 0, we get 



Mb 
Mu 



Bu 



2.5/6 



A/, 



top 



1.5/6 



1/2 1/12 



D "(7 ' 



Bu 



2.5/12 / ^ 1.5/12 



top 



1/4 1/24 



O-D 



The quantity in the square bracket is A ~ 0.203. This gives the value of the bottom mass as to be 4.7 GeV. 
For the strange we have, 



BTr^^logA^,^ 



-Gr 



4.5 



- log - Gu 
at 



4.5 
~6" 



evr^-logA/^p-Gy 



This leads to 



Ms 
Mu 



Bu 



4.5/6 



Mr, 



VAf, 



top 



4.5/6 



-1/2 



Bz 

Bu 



4.5/24 



My 
Af, 



4.5/24 



-1/8 



Here A comes out to be near 0.196 and the strange mass is found to be 0.168 GeV. 
Proceeding further to the down quark, we get 



which gives 



Ma _ (BzV'^'f Mu\ 
Mu \BuJ \MtopJ 



6.5/6 



,-1/2 _ 



Bz 

Bu 



6.5/36 



Mr 



M, 



top 



6.5/36 



-1/12 



A^ 



This gives A ~ 0.19 and Af3=0.0053 GeV. The ratio Mg/Md — 31. However, from eauation (|1.31|) . this ratio 



Ms_ 
Ma 



X-^ = (0.2)-2 = 25. 



These can be considered as gauge interaction corrections to the Wolfenstein parameter. 
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7. THE MASSES OF THE LEPTONS 

Ignoring the non-linear terms proportional to Ap—A, the RG Equation, retaining the top quark only, for the six 
leptons are 



8.^1 log - 


—Ge 




8.^1 log M,! - 


—Ge 




8.^1 log M.^ - 


—Ge 




8.^1 logM.^^ - 


—Gn 


+ 3M2„p 


8.^1 logM.^^, = 


—Gn 


+ 3ACp 


8.^1 logM.^^ - 


—Gn 


+ 3M2„p 



For the first three electron- leptons, we write, with choices based on A ~ .22 and rotational integer np, 
d 



8.^-logA4^ 



8.^-logAf^ 
8^^|logAf,2 



-Ge + I (stt^I logBHr) -Gu)-f^ ^^\ogMl^{t) - Gu 



17 



24 
11 



Ge + ttA^^^. ^ogB'it) -Gu]--:[ St:' - log M'{t) - G 



dt 



17 



-GE + fAS.'j^logB'{t)~Gu 



24 
11 
24 



d 



dt 



[sii'j^\ogMl^{t)-Gu^ , 



and obtain 



Bz 
Bu 



5/4 



5/4 



-1/2 



\ ^/^^ Afi 



Afy 



5/32 

U \ -1/16 
Or 



top 



n 8 



(7.1) 
(7.2) 
(7.3) 
(7.4) 
(7.5) 
(7.6) 

(7.7) 
(7.8) 
(7.9) 

(7.10) 



which gives = 0.209; = 4 x lO""* GeV. Similarly 



M„ 



Mu \Bu J 



By 



17/24 



Mj. 



Mi 



top 



17/24 



-1/2 



which gives A^, = 0.238; Af^ = 0.09 GeV, 
and 



B^ 
Bu 



17/120 / , .r \ 17/120 
/ luij 



top 



-1/10 



- A^ 



(7.11) 



Mr 



BzV'/'U Mu Y'/'' -1/2 



BzV'"' (MuV'"' -1/6 
Bu) \M,,p) 



n 3 



(7.12) 



which gives A^ = 0.23; Mr = 1.36 GeV. 

The masses of the neutrinoes have not yet been measured, only limits have been set. The reported values, as shown 
in Table-HTTI fall into a good pattern, namely, M^,^ = Mu\^^, M^^ = Mu\^ and M^^ = A%A^, A=0.22. So the 
equivalent, degeneracy lifting equations, showing the choices of np are 



8.^-logAf,^^ 



8.^^1ogA.^ 
8.^1 log Af^ 



-Gn + I {&Ml^ -Gu)-2 {dM^^ - Gu) + ^Gu, 



-G 



N 



d 
It 



tt^-^IokB^ ] - 2 [^87r2^1ogAf2 
-G„ + |(8,^|,„gB^)-(8,^|,„gM?„, 



top 



1 



-Gn + ( 8^'^ logS^ ) - ^ ( Stt^-^ log a4p 



d 
di 



'fU: 

Gu, 
Gu- 



(7.13) 
(7.14) 
(7.15) 
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The neutrino masses are obtained as 



N "(7 \ n 



h/2 



Mr, 



1/2 1/ A ( Bz 



\Bi 



3/2 



Mu 



top 



-1/2 1/4 f Bz\( Mu 



1/2 



~ 3 X lO^^GeV, 
~ 1.4 X lO^^GeV, 



^ ~ 1.3 X lO^^GeV. 



(7.16) 
(7.17) 
(7.18) 



This completes our calculation of the masses of the quarks and leptons in terms of known and calculable quantities 
r]F, satisfying the RG Equations. It is absolutely necessary to write down the solutions of the RG Equations for 
each quark and leptons separately, then alone, the contributions from the gauge interactions and the Yukawa mass 
coefficients ApH can be ascertained. 

The values of the coefficients r^p used above have been calculated, upto to the nearest fraction, from the cited 
rotational integers in Table- HllI and the top coupling coefficients Np from the RG Equations. They are given in 
the Table-HVI 



TABLE IV: Coefficients of rip for fermions 



quarks 




Nf 




leptons 


np 


Np 


■riF 


c 


3 


3 





e 


8 





30/24 


u 


7 


3 


1/2 




4 





17/24 


b 


2 


1 


1/4 


T 


3 





11/24 


s 


4 





3/4 




16 


3 


2 


d 


6 





13/12 




9 


3 


1 










Ur 


6 


3 


1/2 



8. RUNNING MASSES OF THE FERMIONS 

The exact one loop solution of the RGE for change in mass values with energy is given by eauation H2.16|) . With 
Zp{t) = Stt^^ log B'^{t), this equation reduces to 



and 



If we take 



and 



ith 



Mlp 



MlpBlitx) , Ap 



77^ / dh 



Ml{t) Ml Blit) • 8^2 ^i(i) ' 

Bp{t)/Bp{tx) 



Mp(t) = Mp{Mx) 



Ml Af ftx Bl(H) 



1/2- 



Bp{t) = A"-(i-*/*-) = exp inp{l - t/tx) log. A) , 



^^^y2np{l-t^/tx) 



1/2 



X 



^^^g-2npti/txlog,A 
tx I 



2np logg A 

Jx l_ 

2np logg A 



g-2n_F log^ A _ ^-2npt/tx loge A 



^— 2nF y^ — ^npt/tx 



(8.1) 



(8.2) 



.3) 



(8.4) 



(8.5) 
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we have, from equation (|8.2|) 

Mpit) ^ Mf{Mx) 



1/2 



(8.6) 



We shall not use this for constructing the tables. A simpler approximate form ignoring Ap terms of equation (|8.6|l . 



■7) 



is presented in tabular form in Table-0to Table- lyTTI This gives the same CKM matrix. Now, we present our results 
in the following tables. 



TABLE V: Variation of Mtop{t) with t= log^iii/Mz) 



t 





3 


6 


9 


12 


15 


18 


21 


24 


27 


30 


33 


Mtop{t){GeV) 


175.88 


166.98 


159.60 


153.25 


147.61 


142.46 


137.64 


133.02 


128. 51 


124.51 


119.56 


115.00 



From the above values, we find that ^^dt'^^^^ is nearly zero. The top mass variation is approximately 



Mtop{t) 



175. - — X 60 
33 



GeV. 



TABLE VL Running Lepton Masses in GeV, tx=33. 



t/tx 


Afe 


A4 


Mr 


A/,, 


A/.„ 


A/., 





0.00051 


0.105 


1.77 


3T0"'-' 


1.9-10-* 


0.018 


0.1 


0.00175 


0.2115 


2.687 


3.43T0"* 


7.19-10"* 


0.043 


0.2 


0.006 


0.426 


4.079 


3.92-10"'^ 


2.72-10"^ 


0.104 


0.3 


0.0206 


0.8575 


6.192 


4.49-10"® 


1.03-10"^ 


0.25 


0.4 


0.02 


1.726 


9.399 


5.1340"^^ 


3.9-10"^ 


0.6 


0.5 


0.242 


3.476 


14.267 


5.87-10-'' 


0.1478 


1.44 


0.6 


0.830 


6.999 


21.658 


6.72-10"^ 


0.56 


3.46 


0.7 


2.849 


14.081 


32.876 


7.68-10"^ 


0.01 


8.3 


0.8 


9.774 


28.369 


49.910 


0.482 


8.02 


19.93 


0.9 


33.527 


57.118 


75.758 


10.05 


30.374 


47.9 


1.0 


115.0 


115.0 


115.0 


115.0 


115.0 


115.0 



TABLE VIL Values of Running quark masses in GeV (other than the top quark), tx~S3. 



t/tx 


Mc 


A/„ 


Mt 


Ms 


Ma 


0.0 


1.27 


0.0042 


4.23 


0.159 


0.0075 


0.1 


2.62 


0.0116 


5.88 


0.307 


0.0196 


0.2 


3.12 


0.0324 


8.18 


0.593 


0.051 


0.3 


4.9 


0.09 


11.39 


1.146 


0.135 


0.4 


7.7 


0.25 


15.85 


2.213 


0.354 


0.5 


12.0 


0.695 


22.05 


4.276 


0.9287 


0.6 


18.96 


1.9307 


30.687 


8.26 


2.434 


0.7 


34.16 


5.3636 


42.696 


15.95 


6.383 


0.8 


46.69 


14.9 


59.406 


30.82 


16.73 


0.9 


73.28 


41.395 


82.651 


59.59 


43.86 


1.0 


115.0 


115.0 


115.0 


115.0 


115.0 
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The graphs are much more reveahng. The variation of mass for the Up quarks, including the top, have been shown 
in figure-^] for the Down in figure-|21 , for the Electrons in figure-|31 and for the Neutrinoes in figure- 2| All have been 
compressed in figure-|Slto allow a glance at the totality of the descent and ascent of the masses to the unification mass 
of 115 GeV. We believe that an exact analysis will not differ much from those presented in these figures. 



9. CONCLUDING REMARKS 

Deo and Maharana^^ have already given exact solutions for the one loop RGEs in MSSM. They proved that all 
the fermions might have originated from a common mass Mu ~ 115 GeV at the GUT energy of Mx — 2.2 x 10^^ 
GeV in a perturbative scheme. This has not been confirmed by authors working in this field as they have not looked 
at the solutions which proved the same wayas one loop gauge sector renormalisation group equation. 

As the energy diminishes, the mass of the top increases to 175 GeV at the Z-meson mass Mz, whereas, masses 
of all other quarks and tiny leptonic masses increase from their value at Mz to Mu of each 115 GeV. Hopefully, the 
exact nature of this variation could be obtained from RGE. By introducing an auxiliary function B[t) through the 
mixing terms of RGE, a very simplified expression can be obtained for both the mass values at Mz{t = 0) and their 
variation till Mx when, the masses of the fermions which keep changing, attain the final value at Mx- 

The success of the analysis of the Wolfenstein and the ratio parametrization of RRR given by eauation H1.31() is 
clearly brought out by our solutions. An extension to running masses at linearity level, leads to a very simple formula 
for the fermions other than the top, 

M,[t)=Mu {^f'^' A-(-^\ (9.1) 

where M|f^' = Mf{Mz) — X"'' ■ The increase is exponential. Much more exact analytic studies are needed to deduce 
the values of the np from RGE. As a first step, we follow up the texture analysis procedure by introducing a similar, 
but not the same function i?(t)'T3'|, which is such that d\ogB{t) — d\ogMtop(t). The resulting non-uniqueness is 
taken advantage of. Introducing known gauge couplings and the values of n^, the effect of gauge interaction have 
been calculated. 

In the foregoing sections 1 to 8, there has been two important omissions. We have not remarked neither about 
higher loop effects nor about the threshhold corrections. They may be very large due to very heavy top. There has 
been no mention of the CKM phases. But, we work with three generations only. So, in this case there is precisely, 
one phase angle </>. This is defined from the Wolfenstein parametrisation as 

1 

2 

(9.2) 

and can be deduced by extending the Gatto-Sarton-Tonio-Oakes (GSTO)^^ to next leading order of solution of RGE, 
A = (Af<i/Ms)i/2. Putting in the numbers np from table-Ill, we get 

cos(/) ~ c± -0.1, </)~95°. (9.3) 

This result is reasonable. 

It may be argued that there are too many parameters still in the guise of the rotation integers np even though A 
could be determined. Let us recall the case of the unification of the strong, the weak, the electromagnetic and the 
gravitational interactions. One consciously omits gravity and lists the other three in order of their strengths. Here the 
top mass, with the largest value, has been deduced from the knowledge of the unification mass and gauge couplings 
quite accurately. For F=2, the two successive rotations i.e., to in gives nF=2 and MjjX'^ falls near the mass of the 
bottom. This process of complete successive rotations yields the masses of all the quarks and electrons, except nF=5. 
However, there is a GUT prediction by Georgi and JarlskogQ, that the ratio Mg/Md ~ 25 ~ A^'^. This may be the 
reason why np jumps from 4 to 6. There is also another accurate prediction 

M,/M, ^ ^ MJM, 




{l-MjMsf ■(1-Me/M^)2- 

With these GUT supplements, there is no other unknown additional parameter left. The unification mass Mjj ~ 115 
GeV is the only input. 
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Corroborating and in consonance with a lot of excellent work in neutrino physics have been done by many workers 
,e.g. Babu, Mohapatra and Barrj23|- The neutrino masses fall with the rotational integers, determined for the 
neutrinoes in the table-V in a very nice way. 

The figure-5 illustrates the significance of this approach of solving the RGE. The similarity of figure-5 with gauge 
unification graph is striking. Thus, we prove, beyond doubt that the original ideas of Pati and Salam^] is correct as 
well as the SUSY standard model. We have given many extrapolatory ideas which can be successfully implemented 
so that SUSY standard model world can be characterised by only three parameters, the GUT mass Mx — 2.2 x 10^^, 
the GUT coupling strength acuT — ^tz common-origin -fermion mass Mu ~ 115 GeV. This is a simplification 
beyond expectation and very much gratifying. 
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FIG. 2: Variation of masses of Down-quarks in Gev with ijtx , t=log(/i/Mz) and tx=33. 
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FIG. 1: Variation of masses of Up-quarks in Gev with t/tx , t=log(///Mz) and tx=33. 
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FIG. 3: Variation of masses of Electrons in Gev with i/tx , t—\og{n/Mz) and tx=33. 




FIG. 4: Variation of masses of Neutrinos in Gev with t/tx , t=log(/u/Mz) and tx=33. 



